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RELATIVISTIC TODA CHAIN AT ROOT OF UNITY III. 
RELATIVISTIC TODA CHAIN HIERARCHY 

S. PAKULIAK AND S. SERGEEV 



Abstract. The hierarchy of the classical nonlinear integrable equations associated with 
^ ' relativistic Toda chain model is considered. It is formulated for the iV-powcrs of the quantum 



operators of the corresponding quantum integrable models. Following the ideas of the 
paper pi[ it is shown how one can obtain such a system from 2D Toda lattice system. 
The reduction procedure is described explicitly. The soliton solutions for the relativistic 
Toda chain are constructed using results of |I2] in terms of the rational tau- functions. The 
vanishing properties of these tau-functions are investigated. 



^ . Introduction 

en 

O ' This paper is a supplement to the papers pl, 0] . It contains proofs of several statements 

O ■ used in these papers and related to the classical relativistic Toda chain (RTCh) model. The 

O \ papers [jl], 0] deal with the quantum integrable model based on the Weyl algebra at the roots 

of unity 

a' 

X. 

H , Namely, the simplest possible model which is the RTCh model at the root of unity was 

considered there. The parameter lo being different from unity describes the deviation of the 
quantum model from the corresponding classical counterpart. 

In the quantum integrable models of this type, in contrast to the models where the quan- 
tization parameter is in general position, there exists a non-standard classical limit. In usual, 
in order to get a classical model from the quantum one, the quantization parameter is sent 
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2 S. PAKULIAK AND S. SERGEEV 

to unity and the Poisson bracket appears in this classical limit as normalized commutator. 
In this limit the quantum observables becomes classical ones, namely, C-valued functions 
depending on the evolution parameters (times or coordinates). In the quantum integrable 
models based on the Weyl algebra at the root of unity there is another possibility. One may 
consider the A^-th powers of the quantum dynamical variables 



(2) 


N N 


Writ 


ing the deformation parameter in the form 


(3) 


^ ^ ^2ni/N_^e/N2 



and considering the limit e — > |^ one may obtain from the commutation relations (|I|) the 
Poisson bracket 

(4) {Um, Wn} = SnmUmWn ■ 

It is clear that logarithms of the dynamical variables form the canonical symplectic form 

(5) n = y^ d \n Um /\ d \n Wm ■ 

The important feature of this unusual classical limit is that the quantum variables do not 
disappear and their A^-th powers being the parameters of the unitary representations of the 
Weyl algebra form the classical integrable system. The quantum analysis of such a combined 
systems becomes more transparent. In particular, in the paper |^ a modified Q-operator 
was constructed for the quantum RTCh model at the roots of unity using several result from 
the corresponding classical integrable system. 

The aim of this paper is to prove some of these results from the point of view of the 
standard theory of integrable nonlinear differential equations. In this approach the integrable 
equation appears usually as a representative of some family of equations sometimes called 
hierarchy of integrable equations. In some cases, these hierarchies are certain reductions of 
a most general system of equations. The subject of the present paper is the RTCh equation 
and the associated hierarchy. It is known that this model can be obtained as a reduction of 



the 2D Toda lattice (2DTL) system |jTT], |T2[ and its integrable discrete analogs 



Relativistic generalizations of the classical integrable equations attracted much attention 



after seminal work by S. Ruijsenaars [in] and have been investigated in several papers (see, for 
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example, 0, |^, |^ ) • The bilinear formalism for one of most simple example of these relativistic 
generalizations, namely, for RTCh equation was presented in [^]. In the latter paper the 
M-soliton solution to this equation were constructed in terms of Casorati determinants. 
We show that these solutions coincide identically with the standard soliton tau-functions 
constructed in KP or 2DTL theories [^, |16|, |T^ providing that parameters of solitons lies 
on the rational curve, which is related to Baxter's curve appeared in the quantum integrable 
models based on the Weyl algebra at the root of unity. 

The paper is organized as follows. In Section |] we consider RTCh model in finite volume 
(with lattice size equal to M) and using L-operator approach construct RTCh hierarchy in 
the framework of slightly modified AKNS construction. In Section S we obtain the RTCh 



hierarchy as certain reduction from 2DTL theory following ideas of the paper |TT[ . In last 
Section the soliton solutions to the first equation of RTCh hierarchy are considered and iden- 
tified with those used in papers ||1|, 0. The Appendix contains a short review of Casoratian 
technique, which helped authors of the paper |]12l to find soliton solutions of RTCh equation. 



1. RTCh hierarchy in t-matrix formalism 

In order to establish relation with the papers [H, 0] we construct in this Section the 
RTCh hierarchy using language of L-operators and r-matrix formalism. We will show that 
RTCh hierarchy can be obtained as slightly modified AKNS system. In particular, we will 
explain the meaning of the constant shift in the corresponding trigonometric classical r- 
matrix mentioned in the paper [||. 

The L-operator for RTCh in the form of 2 x 2 matrix 



(6) Lm 






y/UmWm \ ZWm 



appeared first in |p. In (|^) z is a spectral parameter, k G C is the complex parameter and Um, 
Wmi m E Tj satisfy the local Poisson bracket (^). This L-operatorQ is associated with lattice 
site m. In this section we will consider the lattice of the finite volume, so m = 1, . . . , M. 



^We used here a gauge-equivalent dependence of the L-operator ^ on the spectral parameter in contrast 
to [|l], H and also changed its normalization. 
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It is easy to show p| that the Poisson bracket (H) can be rewritten in r-matrix form 



(7) 



{Lm{z) f Ln{y)} = Snm.[r{z, y) , Ln{z) ® L„Xy)] 



where r{z, y) = r{z, y) + s and 



r{z,y) 



f a{z,y) \ 

(3{z,y) 

(3{z,y) 

V aiz,y) ) 



1 z^ 



y^ 



P{z.y) 



zy 



is the standard trigonometric classical r-matrix for £[2 and 

/ \ 

-1/2 



(9) 












1/2 
/ 



is a constant shift. One may show that after quantization this constant shift is related 
to twisting of XXZ quantum i?-matrix in the commutation relation of the corresponding 
quantum L-operators (see @, 0). In the quantum case the twisting of the quantum R- 
matrix is necessary in order to use minimal finite-dimensional representations of the Weyl 
algebra at the root of unity. On the classical level such a shift of the classical r-matrix 
will simplify the corresponding classical integrable hierarchy by removing higher non-linear 
terms (see discussion around (p^). 

In order to introduce the evolution of L-operator (^) with respect to continuous times one 
have to consider the monodromy matrix 



(10) 



Li L2 



^M 



A{z^) zB{z^) 
zC{z^) D{z^) 



where v4(2^), B{z'^), Clz"^) and D{z'^) are polynomials with respect to z'^ of the orders M, 
M— 1, M—1 and M—2 respectively. The commutation relation (|^ shows that the coefficients 



M 
k=0 



Tfc, k = 0, . . . , M of the trace of the inonodromy matrix ( |I0| ) 

M 

(11) T{z') = A{z') + D{z') = J2Tk 
are in involution with respect to this Poisson bracket 

{Tfc,T„} = 

and can be taken as Hamiltonians which correspond to evolution of the L-operator (^ with 
respect to times t^ 

(12) ^ = {L^,T,}. 



Our goal is to rewrite the equation ([12|) in the zero curvature form. To do this we will 



use a method described in the book ||^ . We have 

(13) {Lra{z),T{y)} = Vm-iiz,y)Lmiz) - L^{z)Vm{z,y) , 
where 

Vm{z,y) = tT2[{id® Li{y))-- -{id® Lm{y))f{z,y) x 

(14) X (id ® Lrr,+i{y)) ■ ■ ■ (id ® LM{y))] 

and tr2{A ^ B) = tr(i?) ■ A. Moreover, since the singular part of the r- matrix (^ is propor- 
tional to the permutation operator the following equality holds 

lim (2:^ - y^)[Vm-i{z,y) Lm{z) - Lm{z)Vm{z,y)] = 

which shows that r.h.s. of (|13|) is polynomial with respect to z and y 

M 

(15) VUz,y) = Y.Qli^\z)y'\ 

k=0 

where 

k 

(16) Q^r^\z) = Y,<^,,{T)Q^^\z), 

p=0 
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, , / ^2" /2 

\ -z2"/2 

(17) V 

^ / /l2n/2 

\ -/l2„/2 ^ 

and the entries of the matrix polynomials Qm (z) h2k, ^2k-i and f2k-i are differential poly- 
nomials of the functions u^ and Wm+i only. The coefficients a^piT) in decomposition (jT^) 
are polynomial functions of the integrals of motions T^. 

We can see that the structure of 2 x 2 matrix polynomials ([17D looks very similar to those 
in the AKNS construction. Let us introduce new time variables xj- such that evolution of 
L-operator (|) is given by the zero curvature equation 

(18) ^^ + Q^ZUz) LUz) - LUz) Q^:i\z) = . 

In (^) we can choose integer n to be any non-negative integer and call the obtained 
hierarchy of equations as RTCh hierarchy. Let us write explicitly several first equations of 



this hierarchy (see the Section y, formula ( piP , where linear combination of the systems ([T 



and (|20|) is identified with RTCh equation): 

- dlnum _ dlnwm _ ^ 

dxQ dxo 



(20) 



OX I 



(21) 



dX2 

d\\iw. 



m 



dxo 



dxi 

= UmWmWm+l{KUm + Mm-l) + 

+ Um^iWm-lWm{K.Um-l + Um~2) + wl^^KUm + Mm-l)^ , 



The compatibility condition between the systems (^) and ( pTj ) falls into quite simple form: 

( 



dUm d'^Urr, „ OUr 



'"m „ ^ ""m 



f22) 



0X2 OXl OX I 

0X2 oxf Ox^ 
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As we already mentioned, the structure of the matrices Qm (0) shows that the whole 
RTCh hierarchy can be obtained as a version of AKNS hierarchy |^, |Tl| with small modi- 
fication related to the shift matrix s ^. In (0) we decompose the matrices Qm into two 
pieces. The first one correspond to the r-matrix (^ and to standard AKNS construction, 
while the second one to the constant shift matrix s (|^). If the parameter e in (^) is free then 
the first equation in the corresponding AKNS construction, written as the zero- curvature 
condition for matrices (p!7|), 

has the form (we renamed Um = e and Wm+i = / for simplicity) 



(24) 



* - S%4„/|i + 2(2.+ l) e^|^ + (.+ l)e3f 



dx2 dxf dx^ \ dx^ 

Sf - «V+4«/|^ + 2(2.H-l)f/'|l-(. + l)eV' 



dx2 dx\ dx^ \ dx^ 



It is seen from (0) that these nonlinear equations are simplified essentially at the choice 
of the parameter e = — 1/2 by removing the higher nonlinear terms. The same phenomena 
occurs with all higher hierarchy equations. We would like to repeat once more that such a 
simplification is a classical analog of the possibility to use the minimal representation of the 
Weyl algebra in the corresponding quantum model. 

Due to locality property {Qm depends only on Um and Wm+i) the AKNS construction 
do not describe discrete time evolution (plsl) . As well as in the usual AKNS hierarchy the 
dependence on the discrete time or the discrete evolution (|l^) can be reconstructed after 
using the Backlund-Schlezinger transformation [l^. The L-operator (|]) plays the role of 
matrix which generates this transformation. It is possible to develop this AKNS construction 
in full details using the methods described in the book 0, but instead we will describe in 
the next section the RTCh hierarchy as a reduction of 2DTL. 

2. RTCh as reduction of 2DTL 
In this section we recall the basic facts on 2DTL model following the fundamental paper 



T^. Then we consider a certain reduction of this system which leads to the RTCh hierarchy. 
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2.1. 2D Toda lattice. In |1J] 2DTL was defined in terms of Z x Z matrices L and L of 
the following form 

(25) L = ^diag[6j(m)]A^' , C = J^ diag[cj(m)]A^' , 

where bi{m) = 1, C-i{m) ^ for any m G Z and Am,k = Sm+i,k is a shift matrix. In 
components, these infinite matrices can be written as follows 

oo 
^m,k = Sm+l,k + / ^b^j{m)Sm-j,k , 

(26) _ ^=° oo 

Cm,k = C-i{m)Sm-l,k + ^Cj{m)Sm+j,k ■ 

For arbitrary infinite matrix A = ^.ggdiag[aj(m)]A'' we define its projections on the 
upper-diagonal (including diagonal) A^ and the strictly lower- diagonal (excluding diagonal) 
A_ matrices by the formulas 

(27) A+ = Y^ diag[aj(m)]A-'' , A_ =Y^ diag[aj(m)]A^' . 

i>o i<o 

2DTL is a system of differential equations for the functions b-j{m), Cj{m), j > and 
c_i(m) with respect to two copies of the continuous parameters (times) x = {xi,X2, . . . } 
and y = {yi, 2/2, • • • } of the form 



f28) 






where Bg = {'C'^)_^_ and C^ = (£ )_. 

The system ( P^ can be obtained as the compatibility conditions of the spectral problems 

(29) CW{x,y) = W{x,y)A, ZW{x,y) =W{x,y)A-^ 

and the linear problems 

dy^Wix,y) = C,W(x,y), dy^W{x,y) = CsW{x,y) , 
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where W(a;, y) and W{x, y) are Baker- Akhiezer functions. Writing them in the form 



^{x,y) = J^diag[i/'j(m;a;,2/)]A ^ ■exp(^(a;,A)) 



i=o 



>V(x,y) = ^diag[^j(m;x,?/)]A^ ■exp(^(y,A ^)) 



j=0 



W{x,y)-' = exp(-e(x,A))-^A--'diag[^;(m + l;x,y)] 



3=0 



oo 



yV{x,y)-^ = exp(-e(i/,A ^))-^A^diag[^*(m + l;a;,i/)], 

i=o 

with iljQ{m;x,y) = ipQ{m;x,y) = 1 and ipQ{m;x,y) ^ for m G Z and using a bihnear 
formahsm one can prove |114[ that the scalar generating series 



-^{m;x,y;X) = A'"exp^(x, A) ^ ^'^(m; x,?/)A ^ , 

j=0 

oo 

^*(m;x,?/;A) = A~™exp^(-x, A) ^ V'*(m; a;, y)A^ 

(31) _ L='_ 

■^{m;x,y;X) = X"" exp^{y, \~^)^'iljj{m;x,y)X^ . 



j=0 

oo 



-^*{m;x,y;X) = X ™exp^(-y, A ^) ^ V^*(m; x,y)A^ 

i=o 

can be expressed in terms of family of tau- function r„(x, y), n E 1^ 

<i/{m;x,y;X) = ^ , 

^*im;x,y;X) = I^A^lI^^ ^ 

(32) >^'''^\^^^ 

Tt7^ \^ 'T^m+i[x,y-e[X)) 
^{m-x,y-X) = ^ , 

"^ {m;x,y;X) = ^ . 

rm{x,y) 

In ( |3lD the function ^(x, A) is 

oo 

e(x,A) = 5^x,A^ 

s=l 

and in ( ^2]) the notation £:(A) means the set {A, A^/2, A^/3, . . . }. 

The formulas ( |5^ and the spectral problems ( POD allow to express all the entries of the 
matrices C and C as certain combinations of tau-function Tn{x,y). In particular, the first 
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coefficient c_i(n) in the decomposition of the matrix L has the form 



(33) 



c_i(m;a;,2/) 



2.2. RTCh hierarchy as reduction from 2DTL. Consider the following substitution for 
the entries of the matrices L and L (j > 0) W% '- 



(34) 



b^j{m) 



Cj[7n) 



c_ilmj 



[K, 



3 3 

j=l j=0 



+ u. 



m—j—l) 1 



"^n^'^+^n^'n^ 



m+iK'^'^m+j ' '^m+j+l) ' 



i=0 



i=0 



""m—1 



Ur. 



In ( |3^ ) we introduced the functions Um{x,y) and Wm{x,y) which depend on the continuous 
times Xs and ys and k G C is a parameter, the same as in the previous Section. 

Comparing the last formula in (^^ with ( p3D we can choose the following parametrization 
of the function Um{x, y) in terms of tau-function 

'Tmi.X.y) 



(35) 



Umix,y) 



r. 



m+l 



'x,y) 



We consider the similar ansatz for the variable Wm{x,y) 



(36) 



Wm{x,y) 



e. 



m+l 



[x,y) 



Om{x,y) 

The formulas (|35|), (|3^) and (|3^) allows to get simple expressions relations for the matrices 
C and £ in terms of RTCh dynamical variables Um and Wm- We have 



Proposition 1. 

(37) 



C 



e 

r 



w 



r 



l + A- 
1 



-U{A-k) 



1-K-^A 
where U , W, T and are diagonal matrices 



U-' (1+A"M W 



T 
6 ' 

-l^ TJ/-1 



T 

e ' 



(38) f/^i = diagf-u 



.±11 



W = diag[u;, 



±11 



T 



±1 ^ diag[r;;^i] , Q^' = diag[C] 



and t/ie rational functions 



1 — X 



in formulas (W^ is always understood as the series 



^ oo 

— — = Vx^' 
1- X ^-^ 
j=0 
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Corollary 1. The matrices C and C which correspond to RTCh reduction satisfy the relation 



(39) 



c-c = c-c 



-K 



Proof. Substituting formulas ( ^4|) into (126|) we obtain 

l-A-'C = (I + k) 
(40) 

'CA + K= = (1 + k) 
where we introduced the diagonal matrices 



e 



r_ 1 + A-i T^ 
'r K- A r ' 



(41) 



r± = diag[r„±i] , 6 = diag[6'„] , 6*^ 



Tm+lTm-l + I'^'^n 



(1 + K)em 

Excluding from formulas (^0[) the diagonal matrix B we can find first the relation (|39|) and 
using this relation in (|40|) we prove the statement of the Proposition. D 

At the first sight the formulas (pTj) does not depend just on variables m„ and Wn, since 
they includes explicitly the diagonal matrices T and 0. But this is not true since, we have 
obvious formulas 



(42) 



e 


A^' 


T 


r 




H 


e 


A^' 


T 


r 




H 



diag 
diag 



n 

.4 = 

n 

i=l 



'^m+i^m+i 



'^m—i'^m—i 



j>0 
j<0 



Using formulas ( pTf ) we can obtain explicit formulas for the elements of the matrices Bs 
and Cs which describe the evolutions of the matrices C and C by the formulas (pH]). This 
will give the evolution equations of the functions Um{x,y) and Wm{x,y) with respect to 
continuous times Xg and ys- Form of these equations is given by the following 



Proposition 2. The substitution (^[j describes the reduction of 2DTL to the couple of the 
RTCh hierarchies with respect to the sets of the times Xg and yg- These hierarchies are given 
by the formulas 

dlnum 



(43) 



dxg 

d In Wr, 
dxg 



E IfiJ.. '"+'-E 



3r 



Jk>0 i=0 
Jk>0 i=0 



r=l 



m — I 



+E> 



r=l 
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and 



r_^s9lnuj 



s-1 



dys 



(44) 

where 
(45) 

and 



dhiw 

\s m 



dys 



i^.>o j=o \ r: 



jfc>o j=0 




Uj(m) = (-)-'if;^ (k(1 - (5j,o)um-i+i + ^im-i) 



(46) 



Wj(m) = (-)%„^ {k{1 - (5j-,o)w-i^.+i + w-i^.) . 

The meaning of the Proposition ^ can be understood as follows. Take the equations from 
the system (|28|) , which describe the evolution of the operator C and C with respect to the 
time Xs- This set contains infinite number of equations, which describe the evolution of the 
entries b_j {m) and Cj {m) with respect to this time. The statement of the Proposition |^ means 
that the whole set of these equations reduces after substitution ( p^D just to pair of equations 
given by the formulas (|^). The same is true for the evolutions with respect to times ys- One 
may verify that first two pair of equations from ( |13| ) which correspond to s = 1,2 coincide 
precisely with equations (|20| ) and (^) obtained from L-operator approach. A general proof 
that the hierarchy (|T8|) obtained in L-operator approach coincide with hierarchy (^) for any 
s should be a repetition of the analogous proof for the equivalence of such a presentations 
in case of standard AKNS hierarchy. We do not go into these details here. 

Explicit formulas (|43|) and (^ij) yield the following 



Corollary 2. The systems (^) and j^^) transforms to each other after the simultaneous 
transformations 

(47) Xs < — > -{-K^ys , Un < — > w~^ . 

The first equation in j^^ ) goes into the second equation of (^^ while the second of 1 ^^ ) to 
first of g^. 
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The transformation (^) allows to forget about dependence of our dynamical variables on 
times ys and consider only the set of equations (^31). 



This symmetry of the reduced 2DTL is related to certain Backlund transformations of 
2DTL investigated in [|T7| and used in the paper []T^ in order to obtain the bilinear formula- 
tion of RTCh equation. Ended, formulas (^0|) contains 'new' tau-function dm given by (l4ll) 



as a certain rational combination of r^ and Om- It is actually unclear from our presenta- 
tion whether the functions dm and 6m are also the tau-f unctions. The last equation of the 
substitution (0) and (|33D allows to identify only r^ with some tau-function of 2DTL. 

Nevertheless, this assertion is true and one can use the following argument to see that the 
function 9m is also a tau-function. It is clear that the symmetry transformation p7|) moves 
the substitution fBll) into another one: 



3 J 

c_i(m) = , 

Wm-l 

where now the ratio of functions Wm should be identified with a ratio of tau-functions ac- 
cording to the formulas (^). The ansatz ( |5BD allows to identify the function 6m with some 
another tau-function of 2DTL. This tau-function should be different since otherwise the 
substitutions (^) and (^Sj) becomes trivial (the product UmWm = 1 for any m G Z in this 
case). 



On the other hand if one uses the substitution (|48|) in order to get expressions for the 



matrices C and C similar to (^OD a new rational combination of tau-functions 



^ ' "^ " (1 + ^)rm 

(k) 

will appear. Let us denote all these tau-functions as a set Tm 

We can extend the value of the index k to be any integer number by repeating the transfor- 
mations given by the formulas ( Pj) and ppj). It was shown in [jT^ that the whole set of these 
tau-functions satisfy a reduced discrete 2DTL (the transformations (ETl) and (|i9|) describe 
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the evolution with respect to discrete parameter k) and all of them are some tau-functions 
of original 2DTL hierarchy. 

(k) 

We will see in the next Section that in the soliton sectors the tau-functions Tm differs 
by the proportional change of the soliton amplitudes (see (p3|)). In terms of the set Tm the 
transformation (147|) can be written in the form 



(50) r!^)^r!^-^)^ keZ 



3. Soliton solutions of the RTCh equation 



m 



In this Section we will give the solitonic solution of the RTCh equation obtained in [12 
terms of rational tau-functions. We will identify them with tau-functions used in the papers 
|jl|, 0] as well as with soliton tau-functions appeared in KP |15, 16| and in 2DTL [14] theories. 



In order to do this we have to transform the RTCh equation (120) into bilinear form analogous 
to one used in the [jT2[. The authors of this paper used so called Casoratian technique to 
find a solitonic solutions to the RTCh equation in form of the Casorati determinants. For 
the readers convenience we will describe part of this technique relevant for our paper in the 
Appendix. 

We start with rewriting the equations (|20|) into bilinear form. Actually, in order to use 
results on the reduction described in the previous section, we have to consider certain linear 



combinations of equations from the RTCh hierarchy. Instead of (pOD we consider a linear 
combination of (0) and (pOD 

■9 In Mm , \ f , ,\ 

— = Wm{K,Um + Um-l) - [H + I) , 

(51) .1 

7{ = Um{HWm + Wm+lj - [K + I) . 

OX I 

such that the system (pT|) has a constant {um = w^ = const) solution. 



Excluding from (^) Wm and replacing Um = e'^"^''^^^^^^^'^^ we obtain (g^ = dqm/dxi 



.52^ " = (Vr + ^ + g.:„_i)(v/r + ^ + gJ _ {VI + k + q'J{Vl + K + q'^^^) 

The latter is the relativistic Toda chain equation. In the non-relativistic limit k — > cxd we 
obtain from (^) the standard Toda chain equation 

„" plm-l—gm. „qm — qm + l 

Hm ^ c 



15 



Proposition 3. |]T2[ The substitution 



'^m\-^lj 



transform the equations ( |31| ^ into bilinear form 

Dl{Tm O Trn) = 2(1 + k) (6',„+i^m_i - T^) , 

(54) -Dl(rm O Tm-l) = (1 + K){TmTm-l - OmOm-l) , 

(i + KjUm.t'm ^ Tm+lTm—l + ^''"m • 

In ( |S^) the Hirota operator D{f o g) = f'g — fg' have been used. 
We can also find a bihnear relation between tau-functions r^ and 6*^ 

(55) -01(6*^-1 O Tm) = (1 + K){0mrm-1 - ^Im-irm) 

which is a consequence of the RTCh equation and Backlund transformation r^ -^ 9^ dis- 
cussed at the end of the previous section. In the Appendix we will show how this bilinear 
equation can be solved using Casoratian technique. 

Let A* and A are two independent complex parameters which belong to the rational curve 

(56) K = A*A-A*-A. 

In the paper [|l| we used the following uniformization of the curve ( ^61) 



A^ = e"^'^ ( a/cos2 (/) + k + cos 

(57) ^ 

A^ = e*"^ f ycos^0 + K + cos 4, 

which is convenient for restricting the model to the chain of the finite size with periodic 
boundary conditions. For real k and the parameters A* and A are complex conjugated. 
Define functions 

(i-A:r(i-A^-i)V^-^p"i 



^ ' (i-a^)'"(i-a;~i)Vi-^^)"i 

f (1 + /t)W2 
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which depend on discrete time m, first continuous time Xi and some additional integer 
parameter k. Define also the M-soliton tau-function by the determinant^ 

^!n^ (xi, {0fc,/fc}|fli) = det ip^^^^_{{xi;(j)jjj) 

\ / «j=i 

One of the statements of the paper [|T2[ is the following 



(59) 



Proposition 4. |T2| The tau-functions which solve the equations ( [5^j are given by the fol- 
lowing specialization of the function^^) 



,M 



(60) 



Tm [Xi; {(t)k, fk}\k=l) = Tm [Xi] {(pk, fk}\tLij , 
Om [Xi; {(j)k, fk}\tLl) = Tm^^ [Xi; {(l)k, fk}\kLl) , 

Om(xi;{^kJk}\k=i) = ^mMxi; {0fc,/fc}|Jlj . 



The tau-functions given by the Proposition || coincide with those used in the papers |jl|, Q 
up to certain normalization factor. Let us change this normalization in order to fix precise 
relation to formulas of these papers and also to demonstrate the connection of these tau- 
functions with standard solitonic tau-functions known from the 2DTL theory [|14|. We obtain 
(A* = A* and A, = A<^J 



(61) 

where 
(62) 

and 
(63) 



-A,r- 



r^' = Il(^: - a;) n Ti^v ^" ^^^-^' ^^L^^ 



2<J 



TL^^ = ^(^*-^*)"'det 



i<j 



A* M 



i,j=^ 



X^^\x,;<j>J) = fe^^*^-^^^^^ 



A* 



m+l 



w)r , <</>) 



a; (1 - A 

A 



1-a;) 



It is clear that factor before Tm in formula (^) cannot spoil the form of the equations (54) 
so can be dropped. Now the tau-functions Tm = 'Tm and 6^ = Tm coincide identically 
with those used in papers [||, If. The only difference is the shift in amplitude (|63|), but this 
is because we used Um = Tm/Tm+i instead of Um = Tm-i/Tm in |l], H- 



^The superscript '^^^ counts the number of Backlund transformations (p5|), not the number of sohtons as 
in the papers (ij ||. 
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The formula ( |62D may be rewritten in the form analogous to the standard presentation of 



solitonic tau-functions appeared in KP ||T3| or 2DTL theories [|T1| with all higher continuous 
times but Xi are freezed: 

M 

A * A 

(64) T^;^^ = det 



6ij - X^^ {xi ; (f)j ,fj)Yl ^hk 



M. ^'-^^ 



*J = 1 



where 



(65) s,,u = s{ct>A'), K^,^0 = |^ ^!_^f 

In order to come from formula ( ]62| ) to (^4|) one have to use explicitly the inverse Vandermonde 
matrix and this is a reason of getting product of factors Sj^k in soliton amplitudes. Actually, 
the determinant presentation (p^ ) is equivalent to the recurrent relation which relates the 
tau-functions with different number of solitons: 

M-l 

fc=i 

where the phase shift is given by cross-ratio 

.... , _ (a:-a-)(a.-a,) 

^ ' " (A*-A,)(A,-A*) • 

Specializing the M-soliton tau-function Tm {xi, {(j)k, fk}\kLi) to the zero value of contin- 
uous time Xi = and setting the amplitudes /^ = 1 we obtain that it becomes proportional 
to 

(67) T^^^O;{^kA}\ti) - det\i-y-' {A:^^--' - A,^^--%^^ ■ 

This proves the following 

Proposition 5. The rational solitonic tau function Tm (0; {0^, 1}|^{^]^) vanishes for values 
of the discrete time m = —M + 1, —M + 2, ... ,0 irrespective the values of the spectral 
parameters (pk- 

This property of the rational tau-function was crucial in [Q in order to construct modified 
Q-operators. 
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4. Discussion 

So far the discrete variable m in the previous sections runs over all integer numbers. It 
is clear how to apply the results of the previous section to finite volume system or how to 
satisfy the periodic boundary condition with respect to discrete time m 

(68) Um[Xn) = Um+MyXji) , Wm[Xn) = Wm+M[Xn) ■ 

One should specialize the spectral variable which uniformize the curve (|56| ) to the finite 
set 

(69) '^" = 1^^ fc = l,...,M-l. 

The case k = or k = M are excluded from (p9| ) because at this value of the parameter 0, 
A = A^ and the solitons disappear. It is clear that in this case there is only finite number 
(M including a zero soliton case) rational tau-functions. 
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Appendix: Casoratian technique 

The goal of this Appendix is to demonstrate Casoratian technique |T^ and prove the 
bilinear relation (^) with tau-function given by the Casorati determinant (|59|). Let us fix 
values of the discrete parameters m and u. The functions ( |55D satisfy the following dispersion 
relations 



r% (m+i/) 

(70) dx, ""^ "'^ ^"'+''+' ' 

V i + « [Vm+k - '^m+k ) - ^m+k+l ' ^ ^ /, . 



im+v) 



19 



Let (^^'' means the column 



(71) 






' (Pm^''\xi] (pi, fi) ^ 



(m+u) I 



\ ^m iXi;(j)M,fM) J 

and consider the following identity for 2M x 2M determinant 



(72) 



^m V^m+M-3 ^m+M-1 Ym+M~2 

^ Vm+M-l Ym+M-2 Ym-1 



—(u-l) 
Vm+M~2 

—iv-l) —[v-l) 

Vm+M-3 V'm+M-2 



The identity is valid since this determinant can be transformed to 



fli:^ ■■■ fm+M-3 





—(I'-l) —(u) 

fm+M-1 Vm+M-2 



—{u~l) 





H^m ^m+M-3 ^m+M-2 



by adding and subtracting the columns and rows. Applying the Laplace expansion to the 
l.h.s. of identity ( |7^ ) we obtain the bilinear relation for M x M determinants 



'rm 'i^m+M-3 ^m+M-l V^m+M-2 



(u-l) (^-l) _ _ (u-l) (u-l) 

Ym~l '■fm 'Pm+M-3 'rm+M-2 



rm Ym+M-3 V^m+Af-2 rm+M-2 



(.-1) (,_1) _ _ (.-1) [u-l) 

rm-1 rm rm+M-3 rm+M-l 



^{u~l) . . .J-~^) ^(-1) ^(-1) 
rm rm+M-3 rm+M-1 rm+M-2 



(.-1) (,_1) _ _ (.-1) H 

rm-1 rm rm+M-3 rm+M-2 



0. 



Using now the dispersion relations ( [70| ) we can identify 



_(!/-!) — (j/-l) —{J^-I) — (I'-l) 

V^m-1 V^m ■ ■ ■ '^m+M-3 Vm+M-2 



—{v-1) — (I'-l) — ('^-l) — (i') 

V^m ■ ■ 'Vm+M-3 fm+M~l fm+M-2 



—(„-!) —(v-l) — (i^-1) —M 

^m ' ' ' '^m+M-3 V^m+M-2 Vm+M-2 



_(i/-l) — (jy-1) — (i^-1) —(1^-1) 

V'm-l Vm ' ' ' '^m+M-3 Vm+M-l 

—{u-l) —{v-1) —{u-l) —{u) 

Vm-l Vm ' ' ' '^m+M-3 Vm+M-2 



T, 



{u-l) 
m—1 ? 



M 



1 dr, 

1 + K dxi 
1 



-H 



VT 



=r. 



{^) 



K 



1 Srl'^-^) 



'm— 1 



\/l + K dxi 
m—1 



SO the bilinear identity for determinants becomes the bilinear identity for tau-functions 



M 



1 dr, 
1 + K dxi 



-H . J'^-'^) 



T. 



m—1 



VT 



=T, 



M 



K 



Vl + fi: t^a;i 



' m m—1 ^ 
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which can be written in the form 

^1 ['m-l " 'm J — K^ ^ i^) y-m ' m-1 ' m-1 'm 

Last equahty coincide with (pSD at z/ = 0. 
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